function u = KernelApply(N, f)
global mvcount;
global Tree;
u = zeros(size(f));
for i = 1:size(f,2)
  mvcount = mvcount + 1;
  u(:,i) = five_solve(N, Tree, f(:,i));
end
